Average formation length in string model 
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The space-time scales of the hadronization process in the framework of string model are inves- 
tigated. It is shown that the average formation lengths of pseudoscalar mesons, produced in semi- 
inclusive deep inelastic scattering (DIS) of leptons on different targets, depend from their electrical 
charges. In particular the average formation lengths of positively charged hadrons are larger than 
of negatively charged ones. This statement is fulfiled for all using scaling functions, for z (the frac- 
tion of the virtual photon energy transferred to the detected hadron) larger than 0.15, for all nuclear 
targets and any value of the Bjorken scaling variable x Bj ■ In all cases, the main mechanism is direct 
production of pseudoscalar mesons. Including in consideration additional mechanism of production 
in result of decay of resonances, leads to decrease of average formation lengths. It is shown that the 
average formation lengths of positively (negatively) charged mesons are slowly rising (decreasing) 
functions of x Bj ■ The obtained results can be important, in particular, for the understanding of the 
hadronization process in nuclear environment. 
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Hadronization is a process that leads from par- 
tons produced in an initial hard interaction to the fi- 
nal hadrons observed experimentally. Two aspects of 
hadronization: (i) the spectra of hadrons produced and 
their kinematical dependences; (ii) the space-time evo- 
lution of the process, at present investigated not equally 
well. First of them has been studied extensively in e + e~ 
annihilation and lepton-nucleon DIS both experimen- 
tally and theoretically. As a result the spectra of had- 
rons and their kinematical dependences are rather well 
known. There are a few successful theoretical models, 
which give transition from initial partons to final had- 
rons through sets of fragmentation functions. The sec- 
ond aspect of hadronization, the space-time evolution 
of the process, despite on its importance, was investi- 
gated relatively little. Although the string model allows 
to study the space-time scale of the hadronization pro- 
cess, only a few works were performed in this direc- 
tion jH 0] . In 10] it was claimed, that for hadrons (as 
composite particles) the very notion of formation length 
is ambiguous because different constituents of a hadron 
can be produced at different lengths. It is then an open 
and model-dependent question which of the two length 
scales plays more important role in hadronization pro- 
cess: (i) constituent formation length l c which is the dis- 
tance between DIS point and the point where first con- 
stituent parton of the final hadron arises; (ii) the yo-yo 
formation length l y which is the distance between DIS 
point and point where two constituent partons of the fi- 
nal hadron meet first time and form particle with quan- 
tum numbers of final hadron but without its "sea". It is 
worth to note that these two length scales connected by 
simple way. 

The string model can give information about space- 
time scale of hadronization, but it gives nothing about 
development of hadronic properties of string during 
hadronization process. At present it is supposed, that 
the unique way to get information about latter is the ex- 



perimental and phenomenological study of hadroniza- 
tion process in atomic nuclei, where string interacts with 
nuclear medium. It is assumed that constituent forma- 
tion length plays more important role in the hadroniza- 
tion process, because beginning with this scale the piece 
of string with a constituent parton on the slow end in- 
teracts with hadronic cross section. As was pointed out 
in Ref. [2] this result follows from the comparison with 
the data on fragmentation of 30GeV pions into p and 
In j4j\ the formation lengths of pions were pre- 
sented in form I = (1 — w)l c + wl y , where w is the prob- 
ability that formation length is l y . Comparison with ex- 
perimental data [5] showed, in case when l c and l y were 
calculated in the framework of standard Lund model, 
that w = 0.1 — 0.17. This result confirmed conclusion 
of Ref. [2] about importance of constituent formation 
length. Further we will consider it as a formation length. 
It is a function of variables v and z (the energy of virtual 
photon and the fraction of this energy carried away by 
the final hadron with energy Eh (z = E^Jv)) and can 
change, as we will see below, in wide region from zero 
to tens (may be even hundreds) femtometers. 

In Refs. [6-8] it was shown that a ratio of multiplici- 
ties for the nucleus and deuterium can be presented in 
the form of a function of a single variable which has the 
physical meaning of the formation length (time) of the 
hadron. This scaling was verified for the case of charged 
pions by HERMES experiment [9]. Now HERMES ex- 
periment prepares two-dimensional analysis of nuclear 
attenuation data using more precise definition of the for- 
mation length of hadron presented in this work. 

In the string model, for the construction of fragmen- 
tation functions, the scalin g fu nction f(z) is introduced 
(see, for instance, Refs. I10h12I1 ). It is defined by the con- 
dition that f(z)dz is the probability that the first hierar- 
chy (rank 1) primary meson carries away the fraction of 
energy z of the initial string. We use three different scal- 
ing functions for calculations: 
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(i) standard Lund scaling function (12 

f(z) = (l + C)(l-zf, 



(1) 



where C is the parameter which controls the steepness 
of the standard Lund fragmentation function; 

(ii) symmetric Lund scaling function il2l,IT7tl 

f(z) = Nz^il - z) a exp(-bm 2 ± /z), (2) 

where a and b are parameters of model, m± = 
-J m\ + p\ is the transverse mass of final hadron, N is 
normalization factor; 

(iii) Field-Feynman scaling function fllll : 



f(z) = 1 -a + 3a(l - z) 2 



(3) 



where a is parameter of model. We will specify parame- 
ters below, when will discuss the details of calculations. 

In the further study we will use the average value of 
the formation length defined as =< l c >. 

The consideration is convenient to begin from L 1 * di- 
rect, Lc ', which takes into account the direct produc- 
tion of hadrons: 



ldWl{L,z,l)/ \ dlD'>(L,z,l) , (4) 



where L = v/k is the full hadronization length, k is the 
string tension (string constant), D^(L, z,l) is the distri- 
bution of the constituent formation length I of hadrons 
carrying fractional energy z. 



D 



h c (L,z,l)=(c^f{z)8{l-L + zL) + 



c h p2 Y,D h cl {L,z,i))e{i)e{L 



i=2 



zL-l) 



(5) 



The functions C pl and C p2 are the probabilities that in 
electroproduction process on proton target the valence 
quark compositions for leading (rank 1) and sublead- 
ing (rank 2) hadrons will be obtained. Similar func- 
tions were obtained in [131] ror more general case of 
nuclear targets. In eq.(5) 8- and ^-functions arise as a 
consequence of energy conservation law. The functions 
D^(L, z,l) are distributions of the constituent formation 
length I of the rank i hadrons carrying fractional energy 
z. For calculation of distribution functions we used re- 
cursion equation from Ref. 0]- 

The simple form of f(z) for standard Lund model 
(eq.(l)) allows to sum the sequence of produced had- 
rons over all ranks (n = oo). The analytic expression for 
the distribution function in this case is: 



D h c {L,z,l) = L{l + C) 



l<- 



(l + zL) 



c+i 



[C h vl 8(l -L + zL) + C^^\ 9{1)9{L -zL-l) . (6) 



In dH (see also Q) eq.(6) was obtained for the special 
case = Cffy = 1. Another special case was consid- 
ered in Ref. Jl5|| . It was supposed, that in electroproduc- 
tion process only the u quarks are knocked out, which 
then turn into observed hadrons. It is clear, that this ap- 
proximation is more or less valid for proton target in re- 
gion of large enough values of Bjorken scaling variable. 
For positively charged mesons it does not differ from 
the first special case, for negatively charged ones it cor- 
responds to special case C pl = 0, C^ 2 = L Nevertheless, 
in this rough approximation it was obtained, that the av- 
erage formation lengths for positively charged hadrons 
are larger than for negatively charged ones. 

Unfortunately, in case of more complicated scaling 
functions presented in eqs.(2) and (3) the analytic sum- 
mation of the sequence of produced hadrons over all 
ranks is impossible. In these cases we limited ourself 
by n = 10 in eq.(5). 

As it is well known, essential contributions in the 
spectra of pions and kaons from the decays of vector 
mesons are expected. But so far, the formation lengths 
of pseudoscalar mesons were considered without taking 
into account this possibility. 

We build the distributions for daughter mesons from 
the ones for parent resonances followed the method 
used in llOl lllll for the construction of fragmentation 
functions for daughter mesons. 

The distribution function of the constituent forma- 
tion length I of the daughter hadron h which arises in 
result of decay of parent resonance R and carries away 

the fractional energy z is denoted D? /h (L,z,l). It can be 
computed from the convolution integral: 



D^ h {L,z,l) = d R ' h 



d 4^ 



,/)x 



f H/H 



(7) 



where 



R/h 



min(l,z/zj in ) and z 



down 



min(l, z/zmax), Zmax (z^^) is maximal (minimal) 
fraction of the energy of parent resonance, which can be 
carried away by the daughter meson. 

Let us consider the two-body isotropic decay of 
resonance R, R — > h\h2, and denote the energy and 
momentum of the daughter hadron h (h = hi or h^), in 

the rest system of resonance, and respectively. 
In the coordinate system where resonance has energy 
and momentum equal and pa, 



y R/h 



R/h 



m R 
1 

m R 



E. 



E. 



(0) 



PR (0) 

E~R Ph 



PR (0) 

E~R Ph 



(8) 



(9) 



where thr is the mass of resonance R. In the laboratory 
(fixed target) system the resonance usually fastly moves, 
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i.e. pr/E r -> 1. 

The constants d R l h can be found from the branching 
ratios in the decay process R — > h. We will present their 
values for interesting for us cases below. 

The distributions f R ^ h (z) are determined from the 
decay process of the resonance R, with momentum p 
into the hadron h with momentum zp. We assume that 
the momentum p is much larger than the masses and the 
transverse momenta involved. 

In analogy with eq.(4) we can write the expression 

for the average value of the formation length L R ^ h for 
the daughter meson h produced in result of decay of the 
parent resonance R in form: 



L R,h 



ldlD R/h (L,z,l)/ / dlD R/h (L,z,l) .(10) 
Jo 



Here it is need to give some explanations. We can for- 
mally consider L R ^ h as the formation length of daugh- 
ter meson h for two reasons: (i) the parent resonance and 
daughter hadron are the hadrons of the same rank (H, 
which have common constituent quark; (ii) according 
to above discussions, beginning from this distance the 
chain consisting from prehadron, resonance and final 
meson h interacts (in nuclear medium) with hadronic 
cross sections. 

The general formula for for the case when a few 
resonances contribute can be written in form: 

L h c = f ldl(a ps D^(L,z,l) + a v y2D R / h (L,z,l))/ 
Jn V r ' 

[ dl(a ps D*(L,z,l) + a v J2D R / h (L,z,l)) , (11) 

J ° R 

where a ps (a v ) is the probability that qq pair turns into 
pseudoscalar (vector) meson. Following Refs. d GJ] 
we use condition a ps = a v = h. 

Let us now discuss the details of model, which are 
necessary for calculations. We will cosider four kinds 
of pseudoscalar mesons 7r + , 7r~, K + and K~ electro- 
produced on proton, neutron and nuclear targets. The 
scaling function Hz) in eq.(l) has single free parameter 
C. It is known [12] that comparison with experimental 
data gives limitation on its possible values C = 0.3 — 0.5. 
Our calculations showed that changing of parameter C 
from the minimal value C = 0.3 to the maximal value 
C = 0.5 leads to the small increasing of the average for- 
mation lengths. Therefore further we will present re- 
sults for standard Lund model obtained for C = 0.3 
only. For symmetric Lund scaling function we use pa- 
rameters G|] a = 0.3, b = 0.58Gey- 2 , and for Field- 
Feynman scaling function 111 ill the value a = 0.77 for the 
single parameter a. Next parameter, which is necessary 
for the calculations in the framework of string model is 
the string tension. It was fixed at a static value deter- 
mined by the Regge trajectory slope fl6l,lr^l 

k = 1/(2ttc4) = lGeV/fm . (12) 



Now let us turn to the functions C pl and C p2 . For 
pseudoscalar mesons they have form.: 



C7T+ 
pi — 



±u(x B j,Q 2 ) + \d{x B j,Q 2 ) 



«1 



pi 



7 9 , 



fu(s B j,Q 2 ) + \d{x B j,Q 2 ) 
T,„=u,d.s e 2 q (q(x Bj ,Q 2 ) + q{x Bjl Q 2 )) 



C7T y"-~»7T 2 

p2 — °p2 — Iq J 



r K+ _ |m(xbj-,Q 2 )7 s + ls{x B3 ,Q 2 )-f q 
pl ~ E q=uAs e 2 M x Bi,Q 2 )+q(x Bj ,Q2))- 



°pl 



^u{x Bj ,Q 2 )j s + \s{x Bo ,Q 2 )^ q 
T, q =u,d,s ^(q(x Bjl Q 2 ) + q(x B j,Q 2 )) 1 



p2 



^p2 



Iq Is 



o 

where x B j = 2 ^ - is the Bjorken's scaling variable; 

Q 2 = —q 2 , where q is the 4-momentum of virtual pho- 
ton; m p is proton mass; q(x B j,Q 2 )(q(x B j,Q 2 )), where 
q = u,d,s are quark (antiquark) distribution functions 
for proton. Easily to see, that functions C pn for hadrons 
of higher rank (n > 2) coincide with ones for second 
rank hadron C pn = C p2 ■ This fact was already used for 
construction of eq.(5). For neutron and nuclear targets 
more general functions C^i = 1, 2) from [13] are used. 
Functions for resonances can be built in analogy with 

the above equations. For example, C P vi = C^ + , where 
i = 1,2. 

All calculations were performed at fixed value of 
v equal lOGeV. Calculations of z - dependence were 
performed at fixed value of Q 2 equal 2.5GeV 2 , which 
give x B j m 0.133. For quark (antiquark) distributions 
in proton the parameterization for leading order par- 
ton distribution functions from [18] was used. We as- 
sume, that new qq pairs are uu with probability j u , dd 
with probability 7^ and ss with probability j s . It is 
followed from isospin symmetry that 7« = 7d = In - 
We use two sets of values for 7: for Lund model Il7ll 
7« : Id '■ 7s = 1 : 1 : 0.3; and for Field-Feynman model 
lu ■ Id ■ 7s = 1 : 1 : 0.5. 

We take into account that part of pseudoscalar me- 
sons can be produced from decays of resonances. As 
possible sources of tt + , ir~, K + and K~ mesons we con- 
sider p+, p°, u>, K*+ and K*°; p~ , p°, w, K*~ and K*°; 
K* + , K*° and <t>; K*~ , K*° and mesons, respectively. 
The contributions of other resonances are neglected. 

In Ref. [10] were presented simple expressions for 
jR/h whjc]-, are c l OS e enough to the experimental data. 
In case of pions we interested in / P//7r , f u ^ and f K I* , 
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on panels c, d; for Field-Feynman model on panels e, f, 
respectively. The values of parameters using in calcula- 
tions are presented also. Of course results of different 
models are quantitatively differ, but qualitatively they 
have the same behavior as functions of z. Therefore fur- 
ther, for illustration, we will use the results of symmetric 
Lund model only. 

Let us briefly discuss why the average formation 
lengths of positively charged hadrons are larger than of 
negatively charged ones. It happens due to the large 
probability to knock out u quark in result of DIS (even in 
case of neutron target). The knocked out quark enter in 
the composition of leading hadron, which has maximal 
formation length. The K + meson has average formation 
length larger than 7r+ meson, because in first case the in- 
fluence of resonances is smaller. The K ~ meson has av- 
erage formation length smaller than 7r~ meson, because 
it is constructed from "sea" quarks and practically can 
not be leading hadron, whereas tt~ meson can be lead- 
ing hadron due to d quark entering in its composition. 

In Fig.2 the average formation lengths for electropro- 



FIG. 1: Average formation lengths for electroproduction of hadrons 
on proton, normalized on L, are presented as a function of z. The 
contributions of direct hadrons as well as of the sum of direct and 
produced from decay of resonances hadrons are presented. Upper 
curves represent formation lengths of positively charged hadrons and 
lower curves of negatively charged ones. Results for standard Lund 
model are presented on panels a, b; for symmetric Lund model on 
panels c, d;for Field-Feynman model on panels e,f, respectively. 



in case of kaons in f K l K and f^/ K . (Sometimes we 
omit the charge symbols in our expressions, but every- 
where it is imply that the different necessary charge 
states for parent and daughter hadrons are taken into ac- 
count.) The function has form /" /7r (z) = 2(1 - z). 
For the other functions common expression f R ^ h (z) — 

-, // R/h R/h\ . i rpi i c R/h j R/h 

l/{z m ax - z m ' m ) is used. The values of z m ' ax and zj in it 
is easily to obtain from eqs.(8) and (9). For instance, for 

z p JZ x « 0.965 



- For K+ and K 

lK— IK- _ 1. j. 

— g, « 



mesons we have 

K' a /K+ _ ^K'°/K- _ 1. 



the p meson decay into pions we receive z„ 
and 0.035. 

For 7T + and it mesons we have d p ^ = |, = 0.3 

and d K 'l* - 1 

d K' + /K+ = j 
d HK „ 1 

In Fig.l the average formation lengths for electro- 
production of different pseudoscalar mesons on pro- 
ton, normalized on L, are presented as a function of 
z. The formation lengths for pions (panels a, c and e) 
and kaons (b, d, f) are presented. The contributions 
of direct hadrons as well as of the sum of direct and 
produced from decay of resonances hadrons are pre- 
sented. Upper curves represent formation lengths of 
positively charged hadrons and lower curves of nega- 
tively charged ones. Results for standard Lund model 
are presented on panels a, b; for symmetric Lund model 



neutron b 
a = 0.3 
\ b = 0.58 GeV : 




FIG. 2: Average formation lengths for electroproduction ofn + , tt~ , 
K + and K~ mesons on different targets in symmetric Lund model, 
normalized on L, as a function of z. 

duction of 7r + , tt~ , K + and K~ mesons on different tar- 
gets in symmetric Lund model, normalized on L, as a 
function of z are presented. It is taken into account, that 
mesons can be produced directly or from decay of res- 
onances. Easily to see, that positively charged hadrons 
have larger L'* than negatively charged ones. This state- 
ment is right even for neutron target. At large enough 
z (for instance z > 0.2) K + meson has maximal average 
formation length and K~ meson minimal. Difference 
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between of positively and negatively charged had- 
rons is maximal on proton target and minimal on neu- 
tron one. It is worth to note, that results for deuteron co- 
incide, in our approach, with results for any nuclei with 
Z = N, where Z (TV) is number of protons (neutrons). 
Average formation lengths of hadrons on krypton nu- 
cleus, which has essential excess of neutrons, do not dif- 
fer considerably from the ones on nuclei with Z = N. 
In Fig.3 the average formation lengths for electro- 
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FIG. 3: Average formation lengths for electroproduction of n + , ir~ , 
K + and K~ mesons on proton target in symmetric Lund model, 
normalized on L, as a function ofxBj- 

production of tt + , tt~ , K + and K~ mesons on proton 
in symmetric Lund model, normalized on L, as a func- 
tion of XBj are presented. It is taken into account, that 
mesons can be produced directly or from decay of reso- 
nances. 

We obtained, for the first time, the average formation 
lengths for different pseudoscalar mesons in the eletro- 
production process on proton, neutron, deuteron and 
krypton. Main conclusions are: (i) positively charged 
pseudoscalar mesons (tt + and K + ) have the formation 
lengths larger than negatively charged ones (it~ and 
K~) on all targets; (ii) contribution from the decay of 



resonances is maximal for tt+ mesons (reach ~ 20% in 
case of symmetric Lund model), for ir~ and K + mesons 
it reaches a few percents, the formation length of K~ 
mesons practically does not feel contribution from res- 
onances; as it was expected, in case of pions maximal 
contribution gives p meson, in case of kaons K* meson. 
It is worth to note that in string model the formation 
length of the leading (rank 1) hadron l c \ = (1 — z)v/ n 
does not depend from type of process, kinds of hadron 
and target. We want to stress, that obtained result de- 
pends from the type of process, kinds of targets and ob- 
served hadrons mainly due to presence of higher rank 
hadrons. Including in consideration hadrons produced 
from decay of resonances diminishes constituent forma- 
tion length. It happens because for producing of hadron 
with fractional energy z we must have resonance with 
z larger than z. The larger is the fractional energy, the 
shorter is the formation length. Of course this statement 
is right for the large enough z (for instance z > 0.2). 

Calculations performed with different scaling func- 
tions: standard Lund B12I1 , Field-Feynman [11] and sym- 
metric Lund lfl2ll showed, that although the numerical 
values of average formation lengths slightly shift, qual- 
itatively they have the same behavior (see Fig.l). 

Which sizes can reach the average formation length? 
At fixed x Bj it is proportional to v. Consequently it will 
rise with v and can reach sizes much larger than nuclear 
sizes at very high energies. Naturally, this conclusion 
is true if the formalism of string model will continue to 
work at such high energies. 

At present the hadronization in nuclear medium is 
widely studied both experimentally and theoretically. It 
is well known, that there is nuclear attenuation of fi- 
nal hadrons. Unfortunately it does not clear, which is 
the true mechanism of such attenuation: (i) final state 
interactions of prehadrons and hadrons in nucleus (ab- 
sorption mechanism); or (ii) gluon bremmstrahlung of 
partons (produced in DIS) in nuclear medium, whereas 
hadronization takes place far beyond nucleus (energy 
loss mechanism). We hope, that results obtained in this 
letter can be useful for the understanding of this prob- 
lem. 
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